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Abstract. Many types of bacteria swim by rotating a bundle of helical filaments also called flagella. Each 
filament is driven by a rotary motor and a very flexible hook transmits the motor torque to the filament. We 
model it by discretizing Kirchhoff 's elastic-rod theory and develop a coarse-grained approach for driving 
the helical filament by a motor torque. A rotating flagellum generates a thrust force, which pushes the cell 
body forward and which increases with the motor torque. We fix the rotating flagellum in space and show 
that it buckles under the thrust force at a critical motor torque. Buckling becomes visible as a supercritical 
Hopf bifurcation in the thrust force. A second buckling transition occurs at an even higher motor torque. 
We attach the flagellum to a spherical cell body and also observe the first buckling transition during 
locomotion. By changing the size of the cell body, we vary the necessary thrust force and thereby obtain 
a characteristic relation between the critical thrust force and motor torque. We present a sophisticated 
analytical model for the buckling transition based on a helical rod which quantitatively reproduces the 
critical force-torque relation. Real values for motor torque, cell body size, and the geometry of the helical 
filament suggest that buckling should occur in single bacterial flagella. We also find that the orientation of 
pulling flagella along the driving torque is not stable and comment on the biological relevance for marine 
bacteria. 
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1 Introduction 

Many bacteria such as Escherichia coli and Salmonella 
typhimurium swim by rotating a bundle of helical flagella 
PP. Nature's simple and ingenious solution for locomotion 
at low Reynolds number has already inspired researchers 
to apply rotating flagella to perform such diverse tasks as 
pumping fluid [2 or manufacturing nanotubes [3]. Even 
artificial helical flagella already exist [4]. 

The flagellum in the bundle consists of three parts; 
the rotary motor, a short and very flexible proximal hook 
that couples the motor to the third part, the long heli- 
cal filament [1,6,5 . The motors are embedded at differ- 
ent locations of the cell wall so that the flagella have to 
bend around the cell body to form a bundle. Bacteria such 
as Escherichia coli and Salmonella typhimurium use this 
bundle to perform a run-and-tumble motion which enables 
them to follow a chemical gradient (chemotaxis). After 
swimming for about Is, the sense of rotation of one mo- 
tor reverses and the attached flagellum leaves the bundle. 
It goes through a sequence of polymorphic conformations 
until the motor reverses its rotational direction again. The 
flagellum returns to its original or normal helical form 
and rejoins the bundle. During this tumbling event, which 
lasts for about 0.1s, the bacterium changes its swimming 
direction randomly. It is interesting that flagella of cells, 
which are clued to a surface, do not bundle [5 ,6 . Fur- 



thermore, each flagellum is reported to be relatively rigid 
and visible deformations due to rotation are not reported 
[5]. In general, videos show a complex behavior of a single 
flagellum when it interacts with other flagella, with the 
wall, or when it goes through different polymorphic con- 
formations [6 j. In this context, recent articles study the 
synchronization and bundling of two or more flagella due 
to hydrodynamic interactions [7,8, 9[fTU] . 

The polymorphism of the flagellum is a fascinating and 

Using a coarse-grained molecular model, we recently ad- 
dressed the question why the normal polymorphic state 
is realized in a flagellum [20] and also developed a sim- 
ple model for flagellar growth [21]. Based on an extended 
Kirchhoff theory for the helical filament, we modeled the 
polymorphism of the flagellum [22] and were able to re- 
produce experimental force-extension curves where a poly- 
morphic transition is induced by an external force [18] . 

In this article we concentrate on the normal form of 
a single bacterial flagellum, model it by the discretized 
version of KirchhofTs elastic-rod theory and develop a 
coarse-grained approach for driving the helical filament 
by a motor torque. A rotating helical flagellum produces 
a thrust force as explained in fig. [I] that adds up along 
the filament and then pushes the cell body forward. We 
report two buckling instabilities of a fixed helical filament 
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Fig. 1. Each segment of a rotating helical flagellum experi- 
ences a frictional force F that is not antiparallel to the local 
velocity v due to the anisotropic friction of a rod. Whereas the 
force component perpendicular to the helix axis averages to 
zero over one helical turn, the parallel component adds up to 
the thrust force. For a detailed treatment see appendix [A] 




Fig. 2. a) An elastic rod buckles under the influence of a com- 
pressional force F and an external torque M. b) The critical 
values F c and M c at which buckling occurs obey a characteris- 
tic relation. The graph depicts relation valid for a rod with 
fixed ends. 



and the applied torque is balanced by frictional forces and 
torques continuously distributed along the filament. Wol- 
gemuth et.al. investigated a rod with one clamped and 
one free end rotating around its axis. They observed two 
regimes separated by a supercritical (i.e. continues) Hopf 
bifurcation. When the rotational frequency exceeds a crit- 
ical value, the straight filament starts to bend and per- 
forms a whirling motion [24] . In Brownian dynamics sim- 
ulations Wada and Netz observed for the same conditions 
a subcritical (i.e. discontinuous) Hopf bifurcation where 
the strongly bent filament nearly folds back on itself [25] . 
On the other hand, a rod tilted with respect to the rota- 
tional axis bends slightly due to friction at low rotational 
velocity. At a critical value, a discontinuous transition to 
a helical rod shape occurs [26] . 

In this article we treat buckling instabilities for the bi- 
ologically relevant helical filament. The problem is more 
complex due to the characteristic rotation-translation cou- 
pling and the fact that we do not fix the orientation of the 
helical filament. The content of the article is organized as 
follows. In sect. [2] we explain how we model the motor- 
driven bacterial flagellum and how we perform the sim- 
ulations. In sect. [3] we present and discuss our numerical 
results for both buckling instabilities for the fixed filament 
and then address the first buckling instability during loco- 
motion of the filament. In sect. [4] we formulate a buckling 
theory for a helical rod and show that it quantitatively 
reproduces the critical force-torque relation in the bio- 
logically relevant regime. We close with a summary and 
conclusions. 



for increasing motor torque. The first instability occurs in 
the biologically relevant regime. The straight helical fil- 
ament starts to bend under the influence of the acting 
thrust force similar to a rod that buckles under its own 
load. The buckling instability is visible as a supercritical 
Hopf bifurcation in the thrust force. It also occurs when 
the filament is allowed to move by attaching it to a load 
particle. We will develop an analytical model based on a 
rigid helical rod that explains the buckling transition and 
reproduces quantitatively the critical force-torque relation 
from our simulations. 

An elastic rod buckles under the influence of a com- 
pressional force and an external torque acting at both ends 
of the rod [fig.^a)]. This is one of the first examples for a 
bifurcation and Euler was the first to provide the theory 
for the critical load force at zero torque. In general, crit- 
ical force and torque for a rod with fixed ends obey the 
relation [23] 



7T 2 =F C 



r 2 



(i) 



where Lq is the rod length and A its bending rigidity [fig. 
[2Vt>)] . Note that eq. (HI) does not depend on the compress- 
ibility or the torsional rigidity of the rod. 

Similar buckling or elastic instabilities occur in the dy- 
namics of rods at low Reynolds number. Here one typically 
applies a torque at one end of the filament. The rod rotates 



2 Modeling the motor-driven bacterial 
flagellum 

We start with a short review of the elasticity model and 
the dynamics of a helical filament in sects. 2.1 and 2.2 
following our previous work [22] and then explain how we 
model the motor-driven hook in sect. 2.3 A summary of 
the simulation parameters follows in sect. |2.4| 



2.1 Elasticity model of a helical filament 

We describe the conformation of the helical filament with 
contour length L by the space curve of its center line r(s), 
where s is the arc length. In addition, we attach a mate- 
rial frame of three orthogonal unit vectors {ei, e2, e%} : to 
each point along the filament so that es points along the 
tangent of r(s) (see fig. [3J. The generalized Frenet-Serret 
equations transport the material frame along the filament, 



d,ei = f2 x e ?; , 



(2) 



where d s means derivative with respect to s. We can char- 
acterize any conformation by the angular strain vector 
17 = j?2j ^3) given in components with respect to 
the local material frame. Along an ideal helical filament, 
spontaneous curvature ^0 and torsion tq are constant. As 
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Fig. 3. The kinematic variables of a slender elastic rod are 
the space curve r(s) of its center line and the material frame 
{ei, e2, 63} attached at each point of the center line. 

material frame for the equilibrium shape of the helical fil- 
ament, we choose the Frenet frame which consists of the 
tangent vector t = e% : the normal vector e\ = n = d s t/Ko, 
and the binormal vector = b = tx n. The strain vector 
then becomes ft = (0, fto,ro). Further parameters of an 
ideal helix are the pitch p = 27rro/(^o + r o) an d radius 
R = ^0/(^0 + r o)- The ratio of pitch and circumference, 
p/2irR — tana, defines the pitch angle a. 

The total elastic free energy of the filament consists of 
two contributions: 




The first term is Kirchhoff 's classical theory for bending 
and twisting, 

/el = ^1? + ^2 " « ) 2 + ^3 " To) 2 , (4) 

where we introduced the bending rigidity A and the tor- 
sional rigidity C [23,27 . Instead of implementing a con- 
straint for the inextensibility of the filament in our sim- 
ulations, we also add a stretching free energy with line 
density 

/st = f(d s r) 2 . (5) 

We choose the spring constant K such that the changes 
in the filament length are below 1.5 %. The filament is 
inextensible to a good approximation. 



2.2 Dynamics of the helical filament 

We mostly performed deterministic simulations, only in a 
few cases we have added thermal fluctuations. We formu- 
late Langevin equations for the location r(s) and intrinsic 
twist (f>(s) of the helical filament. At low Reynolds number 
elastic force per unit length, f el = —dF/dr, and thermal 
force f th are balanced by viscous drag. The same applies 
to the elastic torque per unit length, m el = —ST/Scj) and 



thermal torque m th . Using resistive-force theory, we intro- 
duce local friction coefficients 7y , 7^ and 7^ (see Appendix 
|A| and arrive at the Langevin equations 

[ 7 ||t t + 7± (1 -t®t)]v =f el + f th (6) 
j r uj =m el + rath- (?) 

Here v = d t r is the translational velocity, w = dt4> the 
angular velocity about the local tangent vector t = e%, 
and means tensorial product. The anisotropic friction 
tensor acting on v in Eq. ([6| couples rotation about the 
helical axis to translation and thereby creates the thrust 
force that pushes the bacterium forward as illustrated in 
fig. [l] [28 . Experiments show reasonable agreement with 
the approach of resistive- force theory [29 ,30 j. Finally, the 
thermal force f th and torque ra th are Gaussian stochas- 
tic variables with zero mean, (/ th ) = and (ra t h) = 0. 
Their variances obey the fluctuation-dissipation theorem 
and therefore read 

</th(*, s) f th (t', s f )) = 2k B TS(t - t')5(s - s f ) (8) 
x [711 1 <g)t + 7_l(1 - t <8> t)] , 
(ra th (t, s)m th (t\ s')) = 2k B TS(t - t')5(s - s') lR , (9) 
(m th (t,s)f th (t',s')) = 0. (10) 

In our simulations we use a discretized version of the 
dynamic equations following our earlier work [31,22] (see 
also Ref. [32,33,34 ). We discretize the center line r(s) of 
the filament by introducing 7V+1 beads at locations rW = 
r(s = i • h) and with nearest-neighbor distance h. To every 

bead we attach the material frame {e^\ e!^} [i = 
0, . . . , N) and approximate the tangent vector by 

e 3 - | r(i ) _ r (t-l)|- W 

The transport of the material frame along the filament 
occurs in two steps: First, we rotate about the bond di- 
rection by an angle f2^h to implement intrinsic twist 
plus torsion. Thereafter, we introduce the curvature of the 
filament, by rotating the the bond vector of the ma- 
terial frame about x e^ 1 ^ by an ang 

into the consecutive direction e^ 1 ^. With this procedure 
the free energy densities f cl and / st from Eqs. Q and ([5| 
are discretized and the functional derivatives of the total 
free energy, f el = —ST/Sr and m el = —5F/5(j), reduce to 
conventional derivatives with respect to rW and 0W. In 
addition, we approximate the tangent vector in the friction 
tensor in Eq. g by t«> = (ef + e$ +1) )/\e® + e< <+1) |. 

2.3 The motor-driven hook 

The fiagellum is driven by a rotary motor embedded in 
the cell wall of the bacterium. The motor torque is trans- 
mitted to the helical filament by a short flexible coupling. 
Because of its shape it is called hook. With a well regu- 
lated length of 0.05|Ltm for E.Coli or S.typhimurium and 
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up to 0.1|Ltm for R. sphaeroides it is much shorter than the 
helical filament [35,36,37,38]. It is also shorter than the 
discretization length of h = 0.2 urn which we can employ 
in our simulations as indicated in fig. [4j We, therefore, 
cannot model the hook in full detail. Instead, we repre- 
sent motor and hook by a motor torque that acts directly 
on one end of the filament neglecting the extension of the 
hook. 

Molecular dynamics simulations showed that the hook 
bends and twists easily. This is possible since conforma- 
tional changes of molecular bonds require only a small 
amount of energy [39] . So the hook itself allows the fil- 
ament to nearly assume any orientation outside the cell. 
Hence, it is comparable to a constant- velocity joint. The 
blow-up in fig. [4] illustrates how motor and hook act to- 
gether to drive the filament. The picture also shows the 
rotational degrees of freedom of the filament at the attach- 
ment point to the hook. The filament can rotate about its 
local axis, about the axis parallel to the motor torque, and 
towards or away from this axis. 

The task of the hook is to transmit the motor torque 
to the filament and to guarantee its rotational degrees of 
freedom. In our coarse-grained model, we implement this 
task by balancing all the torques acting on the first mate- 
rial frame {e^°\ e^} that determines the orientation 
of one end of the filament: 

[7^4 0) 4 0) + ^^ 3 (i-4 0) 4 0) )]- (12) 

=M - A^e™ + (f2 2 - *;o)4 0) ] - C(f2 3 - r )ef\ 

The material frame rotates with an angular frequency u>. 
It gives rise to a frictional torque decomposed into a com- 
ponent along the tangent vector and perpendicular to 
it. The length h appears due to the discretization. The fric- 
tional torque is balanced by the motor or external torque 
M = Me 2 , which we assume constant throughout the 
paper, and the elastic torque —5F/8f2. 



2.4 Simulation parameters 

For the bending rigidity we use A = 3.5pNu.m 2 given 
in Ref. [18] as a typical value for bacterial flagella and 
set it equal to the torsional rigidity, C — A. Our previ- 
ous work showed that this is in agreement with experi- 
mental observations [22]. All other parameters are deter- 
mined by the geometry. In our study we use the normal 
state of the bacterial flagellum with spontaneous curvature 
K\ = 1.3/|im and torsion t\ = 2.1/u.m. In the following 
we study a right-handed helical filament although the nor- 
mal state of a real flagellum is left-handed. We calculate 
the local friction coefficients from LighthhTs formulas [40] 
summarized in appendix |a| as 7y = 1.6 • 10 _3 pNs/u-m 2 , 
7_L = 2.8 • 10- 3 pNs/^tm 2 , and 7^ = 1.26 • 10~ 6 pNs, where 
a filament diameter of about 20nm is used. The length of 
the filament is L = 10|Ltm corresponding to approximately 
four helical turns. The discretization length between the 
beads is chosen as h = 0.2 urn. 




Fig. 4. Blow-up: The hook acts as a universal joint between the 
motor embedded in the cell wall and the long helical filament 
which retains its full rotational degrees of freedom. Main pic- 
ture: The hook is much shorter than the discretization length 
indicated by the blue and red segments of the filament. We 
do not model the hook explicitly but let the motor torque act 
directly on the first material frame of the filament which, in 
principle, can assume any orientation in space. 



3 The motor-driven helical filament 



In this section we study in detail the thrust force that the 
motor-driven helical filament generates both when the ac- 
tuated end of the filament is fixed in space or attached to 
a larger load particle, which mimics the cell body. In par- 
ticular, we describe the buckling transitions by illustrating 
the observed filament configurations. 

It is instructive to shortly look at a completely rigid he- 
lical rod first, which does not exhibit translational motion. 
In the low Reynolds number regime, the angular velocity 
uj of the rod and the applied torque M obey the linear 
relation M — Bu>, where B is the rotational friction ten- 
sor. For a long slender helix like the normal form of the 
bacterial filament, one principal axis of B points along the 
helical axis and the eigenvalues in the plane perpendicu- 
lar to this axis are degenerate, in good approximation [40] 
[41]. Now there is a formal analogy to the motion of the 
force and torque less spinning top with axial symmetry in 
classical mechanics [42 ,43 . We just replace the constant 
torque M by the conserved angular momentum and B by 
the moment of inertia tensor. We explain details in ap- 
pendix [B] According to this analogy, the rigid helix in a 
viscous fluids precesses about the constant applied torque 
while also rotating about its helical axis. However, in our 
simulations we observe that as soon as we introduce a fi- 
nite elasticity of the helical rod, the precession is no longer 
stable and the helical filament aligns, for example, parallel 
to the torque. 
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Fig. 5. (a) Thrust force F versus motor torque M. Four different regimes associated with different configurations of the rotating 
filament exist. In regime (iii) and (iv) the minimum and maximum value of the oscillating thrust force are shown. A supercritical 
bifurcation occurs at the critical torque M c \ ~ l.lpNum indicating a buckling transition. A second bifurcation is visible at 
M C 2 ~ 4.2pNum. The red line in regimes (i) and (ii) follows from resistive force theory, (b) Thrust force versus time for specific 
torque values in the four different regimes: (i) M — — l.OpNum, (ii) M — l.OpNum, (iii) M — 1.2pNu-m, and (iv) M — 4.5pNu.m. 
(c) Characteristic snapshots of the helical filament in the four regimes. The red circular arrow and uj indicate a rotation about 
the local helix axis and F l the local thrust force. The green circular arrow and x show the precession about the external torque 
axis. In addition, the trajectory of the tip of the first tangent vector is indicated: (i) The green line belongs to the perpendicular 
orientation of the filament, (iv) red line: fast rotation about helical axis, green line: slow precession about motor torque during 
relaxation of the filament. 



3.1 Force-torque relation and buckling 

3.1.1 Discussion of the basic features 

The motor-driven helical filament creates a thrust force. 
We calculate it as the force component on the first bead 
parallel to the applied torque M = Me z : F = —dT/dr^ • 
We keep here the bead at a fixed position r° — r(0) and 
use the discretized version of the free energv @. Fi smre 
|5|a) plots the resulting thrust force F versus the applied 
torque determined in simulations without thermal noise. 
We discuss the graph in detail. 

A positive torque M produces a thrust force that pushes 
against the anchoring point of the filament. The thrust 
force is constant in time as indicated by the straight line 
(ii) in fig. [5jb). The illustration (ii) of fig. [5^c) shows the 
stable orientation of the helical filament along the torque 
M. It rotates about the helical axis with angular fre- 
quency uj. The local thrust force acting along the helix 
axis is indicated by F l . The tangent of the filament at the 
anchoring point is tilted against M and the tip moves on 
a circle, as indicated by the schematic. Movies for all four 



types of configurations are available in the supplemental 
material. 

A negative torque M generates a negative force that 
pulls at the anchoring point. However, we realized that 
the orientation of the filament along the torque is not 
stable. For long times the filament turns away from the 
torque axis [green arrow in illustration (i) of fig. [5|c)] un- 
e^t-il it reaches a configuration perpendicular to M", where it 
slowly rotates about the local helical axis and slowly pre- 
cesses about M. This motion is also visible for the tip of 
the first tangent vector. The linear increase of F with M in 
the regimes (i) and (ii) in fig. [5|a) fits well with the result 
from resistive force theory for a perfect helical filament, as 
indicated by the line (see appendix [A]). Small deviations 
are visible at higher torques due to elastic deformations 
of the helix which enhance the thrust force. 

At a critical torque M c \ ~ l.lpNumi the thrust force 
starts to oscillate as curve (iii) in fig. [5|c) indicates. Min- 
imum and maximum values of the force are plotted in fig. 
[5|a). They develop continuously from the constant force 
at M c i indicating a supercritical Hopf bifurcation. Illus- 
tration (iii) of fig. [5|c) shows a buckled configuration that 
rotates about the local helix axis with frequency u and 
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precesses with frequency x about the motor torque M 
keeping its shape fixed. The trajectory of the tip of the 
first tangent vector reflects this motion. A straightforward 
explanation is that the helical filament buckles under the 
thrust force generated by the rotating filament. The force 
adds up from the free to the fixed end of the filament 
and puts the filament under compressional tension. This 
is similar to a rod that buckles under its own gravitational 
load [23,27 . In sect. [4] we will develop a theory for this 
buckling transition which is quite involved. Finally, at a 
critical torque value of M C 2 ~ 4.2pN|Ltm a second bifur- 
cation occurs in the force-torque relation of fig. [5^a) . The 
buckled state itself becomes unstable, visible by the fast 
oscillations of the thrust force in fig. [5^b). The buckled 
configuration is compressed until the fixed end becomes 
perpendicular to the motor torque. At this point the fast 
rotation about the local helical axis stops and the thrust 
force averaged over one fast period is approximately zero. 
Now the strongly bent configuration of the filament re- 
laxes slowly and precesses about the applied torque M 
[second configuration in fig. [5jb)(iv)]. The thrust force on 
the anchoring point slowly increases. When the filament is 
sufficiently relaxed, it starts again its fast rotations about 
the local helix axis and the whole cycle repeats. 



3.1.2 Discussion of additional features 

The reported supercritical Hopf bifurcation is also visible 
in other quantities besides the thrust force. We discuss 
here additional properties of the motor-driven helical fil- 
ament. 

To quantify the stability of the filament aligned par- 
allel to the motor torque axis, we recorded the temporal 
evolution of the elastic free energy starting from a small 
disturbance of the aligned state and fit it to the form 



| J 7 - J* 1 « 5F° exp(At) sin(c^). 



(13) 



Here uj is the angular velocity of the rotating helix leading 
to oscillations in T and A is the reorientation rate. The 
result for A is plotted in fig.^a). For positive M below the 
critical torque, the negative A indicates the stable aligned 
state. For small M a reorientation of the filament could not 
be detected within the simulation time. Frictional forces 
are small and hardly deform the helix which, therefore, 
just precesses about the applied torque. Nevertheless, to 
record the thrust force-torque relation, we always started 
from an aligned state at M = lpN|Ltm and then changed 
the driving torque to the desired value and let the elastic 
free energy relax to its stationary value, where we finally 
recorded the thrust force. The small positive A for M < 
indicates the slow reorientation of the filament towards the 
perpendicular configuration The supercritical Hopf bifur- 
cation is located where A changes sign from negative to 
positive. 

Figure [6jb) shows the angular frequency uj for rota- 
tions about the local helix axis as a function of M. The 
linear regime belongs to the aligned state, deviations from 
it occur in the buckled state. The precession frequency \ 
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Fig. 6. (a) Relaxation rate A of the elastic free energy ver- 
sus applied torque M for a small disturbance of the aligned 
state where the filament is parallel to the torque direction, (b) 
Angular velocity uj and precession frequency \ versus torque 
M. The supercritical bifurcation at M c \ is clearly visible. The 
red line is calculated with resistive force theory, (c) Mean end- 
to-end distance (r) in units of the helix length Lo versus M. 
The red dots are results from Brownian dynamics simulations 
with thermal noise included, (d) Standard deviation cr(r) of the 
end-to-end distance in units of ao — R/V2 versus M. Thermal 
noise (red dots) leads to fluctuations about the mean value. 

for rotations of the whole filament about the torque axis 
is plotted in the inset. A non-zero x corresponds to the 
buckled state. 

Finally, figs. [6jc) and (d) plot the mean end-to-end 
distance (r) of the helix and its standard deviation a as a 
function of M, respectively. They are defined as 



(r) 



lim 



(r)) 2 ). 



L)-r(s = 0)\t 



(14) 
(15) 



Whereas (r) is continuous at both bifurcations, the stan- 
dard deviation displays a pronounced discontinuity at the 
second bifurcation in agreement with the behavior of the 
thrust force. We write a in units of <to = R/V2, where 
R is the helix radius, <jq is the maximum value of a in 
regime (iii) where the buckled helix has a constant shape 
but the free end of the filament rotates on a circle with 
radius R. The strong increase of a in regime (iv) is due to 
the oscillating buckled state. 

The rotating filament also experiences thermal forces 
due to the viscous environment. However, since the persis- 
tence length A/ksT w 1mm calculated from the bending 
rigidity A is much larger than the filament length of 10/im, 
we do expect that our results are robust against thermal 



Reinhard Vogel, Holger Stark: Motor-Driven Bacterial Flagella and Buckling Instabilities 



7 



fluctuations. This is confirmed by the end-to-end distance 
(r) in fig. |6Fc) (red dots) which agrees with the determin- 
istic simulations. The standard deviation in[6jd) indicates 
some fluctuations. Below the buckling transition we can 
directly connect them to compressional fluctuations us- 
ing the spring constant of the helical filament, A/(R 2 L), 
calculated in our earlier article [22 . The equipartition the- 
orem gives a /do ~ 0.15 in good agreement with the sim- 
ulated value of 0.18. In the buckled state, the helical fil- 
ament has more opportunities to fluctuate around which 
explains the further increase of a. Furthermore, we observe 
strong fluctuations of the thrust force in our simulations 
which result from the delta correlated stochastic forces 
acting on the fixed first bead of the discretized filament. 
An average over these fluctuations agrees with the deter- 
ministic case (data not shown). The fluctuations will also 
be smoothed out in an experiment which performs some 
temporal average during measurement. 

3.2 Buckling instability during locomotion 

So far we have studied the situation where one end of the 
filament is fixed in space so that it cannot translate. How- 
ever, rotating flagella push the cell body of a bacterium 
forward so that it moves. We mimic this scenario by at- 
taching the filament to a bead of radius a which, for sim- 
plicity, can only move along the z direction. The thrust 
force F generated at the attached end of the filament 
is then used to push the sphere forward acting against 
the Stokes friction force. We observe similar thrust force- 
torque relations as for the case of a fixed filament. The 
aligned state is again unstable for negative torque and 
possesses a larger reorientation rate which might have bi- 
ological relevance as we discuss in sect. [5] For positive 
torque, the aligned state is stable and the thrust force 
grows linearly in the driving torque M until the Hopf 
bifurcation occurs at a critical value M c \ indicating the 
buckling instability. 

Figure [7^a) shows the critical torque M cl as a func- 
tion of the inverse bead radius 1/a. From 1/a = 0, which 
corresponds to the fixed filament, the critical torque in- 
creases linearly in 1/a and then at a -1 « 5/2 turns into 
a slow growth towards the value for the freely swimming 
helix, i.e., 1/a — » oo. In the biological relevant case with 
the cell body size a « 0.5 • • • 2|im, the linear dependence 
of the critical torque on 1/a can be derived based on the 
fact that the critical thrust force F& is nearly constant, 
as we show in fig. [7^b). So the velocity v = F c i/(67rr]a) is 
so slow that the buckling transition is hardly influenced 
by the motion of the helical filament with the attached 
bead. Now, force and torque on the helix depend linearly 
on velocity v and angular velocity uj (see appendix [A| . 
Eliminating uo and setting v = F c i/(67rr]a) at the buckling 
transition, one arrives at 



Cm 



A \\ B \\ 
Cii 



6ttt] a 



(16) 



Here An, Bn, and Cn are the translational, the rotational, 
and coupling friction coefficients parallel to the helical 
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Fig. 7. Buckling transition for a helical filament attached to 
a bead of radius a that can move along the z direction, (a) 
Critical torque M c \ as a function of inverse bead radius 1/a. 
Inset: Blow-up for the biologically relevant regime, (b) Critical 
force F c i versus critical torque M c \. 



axis, respectively. This formula with the coefficients calcu- 
lated by resistive force theory (see appendix|A| reproduces 
the linear increase for small 1/a, as demonstrated by the 
red line in the inset of fig. (7^a). 

In fig. [7^b) we plot the critical thrust force versus the 
critical torque. For biologically relevant values M c \ be- 
tween 1 and 2pN|Ltm the critical force is indeed nearly 
constant. It only shows a very slow linear increase since 
frictional forces due to the motion of the helix stabilize it 
against buckling. At M c \ & 4pN|Ltm the behavior changes 
dramatically. The critical thrust force goes to zero propor- 
tional to M^y (see dotted line) following the behavior of a 
rod that buckles under an applied force and torque as de- 
scribed in the introduction. In this regime the supercritical 
Hopf bifurcation becomes subcritical and hysteresis oc- 
curs. So whereas for small torques buckling is hindered by 
locomotion, for large torque the typical quadratic depen- 
dence F c i M^y is observed. In the following section, we 
develop a theory to describe the observed buckling tran- 
sition. 



4 Buckling theory for a helical rod 

The goal in this section is to formulate a theory that re- 
produces the force-torque relation in fig. (7^b) for the first 
buckling transition of the helical filament as obtained by 
our simulations. Clearly, this relation cannot directly be 
explained by the theory of a thin elastic rod that buckles 
under the influence of an external force and torque which 
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M 



Fig. 8. The helical filament is approximated by a thin helical 
rod of length Lo that is characterized by the effective bending 
rigidity A e R and local friction coefficients of the helical fila- 
ment. The applied torque Me z generates the frictional torque 
m and forces f\\ and f±. 

we shortly mentioned in the introduction in Eq. 0. There 
are several reasons for this. First, the helical filament is 
not just a simple elastic rod. Second, the external force 
that puts the helix under tension is generated locally by 
the rotation-translation coupling of the helix and accu- 
mulated along the filament similar to a rod that buckles 
under its own gravitational weight. Third, the whole fil- 
ament moves with a constant velocity which leads to ad- 
ditional frictional forces and it also precesses about the 
external torque in the buckled state. In the following we 
formulate a model based on the theory of a thin elastic 
helical rod, derive from it a force-torque relation for the 
buckling transition, and compare it to fig. [7^b). 



4.1 Model equations 

To set up our model equations, we approximate the helical 
filament by a thin helical rod where the helicity comes in 
through the rotation-translation coupling in the friction 
matrix. The length of the rod, Lo = sinaL, agrees with 
the height of the helix, where a is the pitch angle. In 
engineering science the buckling of helical springs is a well 
known problem. If the height of the spring is larger than 
its radius, one approximates the spring by a soft rod with 
effective bending, shear, and compressional rigidity [44[ 
1451 146]. In our case, in contrast to classical helical spring 
theory, the pitch of the helix is much larger than its radius. 
We therefore had to generalize the theory of helical springs 
in Ref. [46] to derive an effective bending rigidity of the 
helical rod in terms of the bending and torsional rigidity 
of the filament: 



1 111/ 2 A o 

- — = - — 1 + sin a + — cos a 

A eff 2 A sin a \ C 



(17) 



Details of the derivation are given in appendix [Cj 

To address buckling of the helical rod, we start with 
the balance equations for force and momentum acting on 
a thin elastic rod [23,27 and neglect any inertial contri- 
bution in the low Reynolds number regime: 



F' + f = 0, 
M 1 + e 3 x F + m = 0, 



(18a) 
(18b) 



where ' means derivative with respect to the arc length 
s and es is the local tangent vector. Here F and M 
are internal elastic forces and torques acting along the 
rod, whereas / and m denote, respectively, external force 



and torque densities due to the applied motor torque and 
friction with the surrounding fluid. In addition, bound- 
ary conditions are necessary. At the free end of the rod 
(s = Lo) no external force and torque act, so elastic force 
and torque have to vanish. The end attached to the sphere 
can only move in z direction with velocity v. The exter- 
nal torque Me z is balanced by the elastic torque M(0) 
and the thrust force on the sphere F = Qirrjav equals the 
elastic force at the leading end of the rod, F z (0): 



F z (0) = F = 67rr]ave z 
F(L ) = 0, 



M(0) = Me z 
M(L ) = 0. 



(19a) 
(19b) 



After setting up the problem, we have to explain how 
the different forces and torques entering Eqs. (18) look 



like for the helical rod close to the buckling transition. 
The elastic torque M is proportional to the angular strain 
vector f2 written in components with respect to the local 
material frame {ei, 62,63}: 



M = A eff i?iei + A eff i? 2 e2 + C e ^Q^e^ 



(20) 



where A e fi is the effective bending rigidity of Eq. ( |17| ). 
Since buckling theory considers local displacements ofthe 
rod only, the torsional term and the actual value of the 
effective torsional rigidity C e ff are not important. The for- 
mulation for M is in full analogy to our presentation in 
sect. |2.1[ only the spontaneous curvature and torsion are 
zero for the helical rod which serves as an effective repre- 
sentation of the helical filament. In setting up linearized 
equations in the vicinity of the buckling transition, the 
elastic force F is only needed for the unbuckled straight 
rod oriented along e z . Since the external force density / 
is constant for t he st raight rod, as we argue in the next 
paragraph, Eq. ( |18a[ ) and boundary conditions (19) give 
the linear force profile 



F(z) = /||(L - z)e z with /|| = F/L , 



(21) 



where we introduce f\\ as thrust force F divided by the rod 
length Lq. We will use it as one parameter in the following. 

The straight filament moves with a constant veloc- 
ity ve z and rotates with a constant angular velocity ue z . 
They result, respectively, in a constant frictional force den- 
sity f\\e z and a torque density me z with 



f\\ = a\\v + c\\uj, 
m — c\\v + b\\U). 



(22a) 
(22b) 



where the frictional coefficient cy couples translation to 
rotation. Appendix |A| gives the coefficients ay, by, and cy 
for the helical rod inTerms of the parameters of the helical 
filament. In Eq. ( j2l| ) we have already linked /y to the 
thrust force F. From Eq. (18b) and boundary conditions 
(19), we also deduce a linear torque profile 



M(z) = m(L - z)e z with m M/L , (23) 

where we relate m to the applied motor torque M divided 
by the rod length Lo- So, m is the second parameter in 
our problem. 
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The buckled rod after the first buckling transition in 
our simulations has a constant shape. It rotates about the 
local tangent vector with angular velocity uoe 3 and pre- 
cesses with angular velocity x about the axis of the applied 
torque leading to a local velocity x e z x r. Furthermore, the 
filament translates with velocity v along the z-direction 
and the total local velocity amounts to v = ve z J rX e z x r - 
In the vicinity of the buckling transition, deformations are 
small and in leading order we can identify v and uj with 
the values of the straight rod. Then, the frictional torque 
along the local tangent vector is 



m = me 3 , 



(24) 



where m is already given in Eq. (22b). So, close to the 
buckling instability we can identify m with the applied 
motor torque as in Eq. (23). The frictional force density 
becomes 

/ = /||P||e z + f±P±e z + a ±X P±(e z x r), (25) 
where we use the projectors 



Pii = e 3 (g) e 3 and P± = 1 - e 3 ® e 3 



(26) 



on the directions parallel and perpendicular to the tangent 
vect or e 3 . The force density fn has already been given in 
Eq. ( [22a] and 

f± = a ± v (27) 

characterizes the frictional force density generated perpen- 
dicular to the local rod axis when the rod moves with ve- 
locity v. Since the frictional coefficient a± is larger than an , 
f± acts against buckling. Finally, a±x is the friction due to 
the precession of the rod. We note that a term Py • (e z x r) 
does not appear since it does not contribute in leading or- 
der to f\\ . We also did not include the rotation-translation 
coupling perpendicular to e 3 since the two terms cancel 
each other in the equations, we formulate in the following. 

We will analyze the buckling transition by first con- 
sidering the four parameters m, f\\, f±, and \ as indepen- 
dent and then apply our results to reproduce the force- 
torque relation of the helical rod. Buckling occurs when 
the straight solution r(z) = (0,0, z) of Eqs. (18) becomes 
unstable and a new non-trivial solution occurs at a cer- 
tain parameter set. We, therefore, use the ansatz r(z) = 
(X(z),Y(z), z) and seek two equations linear in X, Y", 
and its d erivatives. We start by taking the derivative of 
Eq. (18b) and use F' = — / to arrive at 



M" 



F-e 3 xf 



0, 



(28) 



where we insert the concrete formulas for Ai", F, /, m. 
We linearize these resulting equations using the identi- 
ties J?i « -y", J? 2 « X", P v e z = e 3 « (X',y',l), 
P ± e z ^ -(X',y',0), and P ± • (e z x r) « (-7,1,0), 
and ultimately arrive at 

= -Y"" + d z (m(l-z)X") 

-f\\(l-&)Y"-f ± Y' + xX, (29a) 
=X"" + d z (m(l - z)Y") 

+ f {l (l - z)X" + f ± X' + X Y. (29b) 



Here we introduced the rescaled coordinate z = z/Lq 
and the dimensionless parameters m = mL\jA Q ^ /y = 

/yLg/Aeff, f± = / ± Lg/A eff , and x = X«±^oMeff- 

Equations ( [29] ) are quite general and several related 
problems follow from them. When forces f\\ and f± van- 
ish, they describe the writhing instability of rotating rods 
[24] . For zero torque and precession frequency, and fn = 

— f± = f z , one arrives at the classical example of a col- 
umn that buckles under its own weight [23][27]. A similar 
problem occurs for microtubuli that buckle under the ac- 
tion of molecular motors [47]. In our case, the force den- 
sity f\\ that causes buckling points along the rod axis and 
f± stabilizes the straight rod for non-zero v. In compari- 
son, the column under gravity always experiences a force 
density along the vertical which gives a force component 
perpendicular to the rod as soon as it buckles and thereby 
supports buckling. 

We complete the linearized dynamic equations ( [29] ) by 
writing the boundary conditions (19) in linearized and re- 
duced form: 



X(0) 


= 


r(0) 


= 0, 


(30a) 


X"(0) 


= -mY'(0) 


r"(o) 


= ml'(0), 


(30b) 


X"(l) 


= 


r"(i) 


= 0, 


(30c) 


X"'(l) 


= 


Y"'(l) 


= 0. 


(30d) 



The first line means that the attached end of the rod can 
only move in z direction and not along the x and y axis. 
The second line means that a torque does not act per- 
pendicular to the z-axis. So if the rod starts to buckle, 
the local torque me 3 has to be equilibrated by a bending 
moment. The free end of the rod is torque less and, there- 
fore, the rod does not bend, as expressed by the third line. 
Finally, the free end is also force free and the fourth line 
follows from Eq. ( 18b ) by setting F = 0. 

To search for nontrivial solutions of Eqs. (J29| in our 
parameter space and thereby identify the buckling transi- 
tion, we procee ded as foll ows. In addition, to the bound- 
ary conditions (30a) and (30b), nontrivial solutions of the 
buckling equations p9] ) can be characterized by X r (0), 
y'(0), X'"(0) and y^O). The principal idea is to use 
them to generate solu tions of Eqs^ J29l) and to fulfill the 
boundary conditions 



(|30cJ) and ( |30df at the free end by 
varying them. However, since X (0) and Y'(0) just deter- 
mine the amplitude of a bent configuration and merely fix 
the rotational degree of freedom about the z axis, they can 
be chosen arbitrary. Instead, we vary two of our four pa- 
rameters, f± and x> to fulfill the four boundary conditions 
at z — 1. As a result, for given rh and f\\ , we determine pa- 
rameters f± and x f° r which non-trivial solutions of the 
buckling equations exist and thereby identify the mani- 
fold of bifurcation points in our four-parameter space. We 
discuss it in the following section. 

4.2 Discussion 

Figure |9|a) plots the manifold of bifurcation points. To 
each parameter triple (m, /y, f±) belongs a specific value 
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Fig. 9. (a) Manifold of bifurcation points in the parameter space (m, f\\, f±). To each parameter triple belongs a specific value 
of the precession frequency %. (b) Buckling curves /y (ra) for different values of the perpendicular force ranging from f± = in 
steps of 25 to 300. The red dots are the critical forces and torques from fig. [TJb). (c) The plane pictures relation (31) between 
/||, /_]_, and m for the helical rod with constant friction coefficients. Intersecting it with the manifold of bifurcation points 
gives the buckling curves in (d). Full blue line: for constant friction coefficients of the undistorted helix, dashed green line: 
torque-dependent friction coefficients. 



of the precession frequency x which we do not discuss 
further here. At positive f± and for small m and h the 
straight configuration of the helical rod is stable. If we 
change the sign of f± , a bifurcation occurs which we inter- 
pret as an instability of the straight rod when it reorients 
towards the perpendicular configuration. We saw this in- 
stability in our simulations when we reversed the driving 
torque as discussed in sects. 3.1.1 and 3.2 Here we keep 
the direction of the torque but reverse the sign of the ve- 
locity v and thereby the sign of f± in Eq. (27) by reversing 



the chirality of the rod. The main result is the surface in 
dark yellow that belongs to the first buckling transition 
observed in our simulations, so at large m the rod is buck- 
led. Finally, at f± ~ and large rh a transition between 
two different configurations of the buckled rod occurs. An 
interesting feature is the ridge in the bifurcation surface. 
However, we could not determine any dramatic changes 
in the buckling of the helical rod close to this ridge. 



Figure |9|b) shows buckling curves f\\(rh) for different 
values of the perpendicular force ranging from f± = 



in steps of 25 to 300. At f± = the typical parabolic 
curve of Eq. occurs. At constant but small value of 
m, the critical force f\\ increases strongly with increasing 

f±. Likewise, one needs large forces f± to stabilize the 
straight helical rod at high torques. The red dots are the 
critical forces and torques from fig. [T^b) determined in our 
simulations. We plot them in reduced units where we cal- 
culate A e ff from Eq. (17). Note that the buckling curves 
develop a shoulder at m around 15 for increasing f± due 
to the ridge in the manifold of bifurcation points in fig. 
[9ja). The two simulation points at large m are close to 
this ridge. We speculate that the transition from a super- 
critical bifurcation observed in our simulations at low m 
to a subcritical bifurcation at large m is connected to the 
existence of this ridge. 



In the rotating helical filament or helical rod, the forces 
/|| and f± and the torque m are related to each other by 
Eqs. (22) and (27). Eliminating velocity v and angular 
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frequency cj, one arrives at 
a± 



aubu 



(hf\\ 



cum 



(31) 



We give the friction coefficients in terms of the helical pa- 
rameters in appendix [Aj Relation (31) defines a plane in 
the parameter space (m, /n, f±) which we intersect in fig. 
|9|c) with the manifold of bifurcation points. The result- 
ing bifurcation curve is then plotted in fig. (9^d) as full 
blue line. At low m we have a remarkable quantitative 
agreement with our simulations (red dots) but we miss 
the slight increase of the critical force fu . We already dis- 
cussed that the location of the buckling transition is sen- 
sitive to small variations in the parameters. We also men- 
tioned in the discussion of the simulation results that close 
to the buckling transition the helical filament is slightly 
deformed. In Ref. [22 we calculated the effective spring 
constant A/(R 2 L) for the helical filament. It gives a rel- 
ative compression of the filament of about 2% for critical 
forces of lpN at the buckling transition, which is negligi- 
ble. On the other hand, we apply a torque along the he- 
lical axis. As a result, one end of the helix twists against 
the other end by an angle A(j)/L = M/A, where we set 
A = C [23]. One end of the helical filament is free, so 
the twist increases linearly from the free end to the other 
and the average value is Mj {2 A). Due to the twisting, the 
radius of the helix changes. One can show that for the av- 
erage twist angle M/(2A) the inverse radius R~ x changes 
to Rr 1 + (cos a) -1 Mj (2 A). Here we keep the pitch an- 
gle a constant, which is confirmed by our simulations. 
The helical radius directly influences the frictio n coe ffi- 
cients bu and cm in Eq. (31) [see Eqs. (38) and (39b) in 
appendix [A] and Eq. (31 ) oecomes a nonnnear function in 
m = M/Lq. Intersecting it with the bifurcation manifold 
gives the green curve in fig. |9jd) which nicely reproduces 
the critical force-torque relation for m < 10. Our theory 
also gives the strong decrease of the critical force fu at 
large m. However, in the effective model the bifurcation is 
shifted to larger torque values. This might be related to 
the ridge in the manifold of bifurcation points. Neverthe- 
less, considering the fact that we approximate the helical 
filament by a rigid rod whose helicity comes in through 
the friction coefficients, we obtain a very good agreement 
with our simulations. 



5 Summary and conclusions 

Bacteria move forward by rotating a bundle of helical flag- 
ella which creates a thrust force that pushes against the 
cell body. In this article we have modeled a single flagel- 
lum based on the discretized version of Kirchhoff 's elastic- 
rod theory and developed a coarse-grained approach for 
driving the helical filament by a motor torque. When in- 
creasing the motor torque, the thrust force reveals a su- 
percritical Hopf bifurcation due to buckling of the helical 
filament. When the torque is further increased, a second 
buckling instability occurs. The Hopf bifurcation is also 



visible when we attach the flagellum to a spherical parti- 
cle, which mimics the cell body, so that the whole model 
bacterium moves forward. Via the size of the cell body we 
can tune the thrust force pushing against the cell body 
and the critical torque for buckling changes. This results 
in a characteristic diagram critical force versus torque for 
the buckling transition [fig. |7jb)]. 

We have developed a theory for the observed buck- 
ling transition by approximating the helical filament by a 
helical rod with an effective bending rigidity and the char- 
acteristic rotation-translation coupling. The basic picture 
is that the filament buckles under the frictional forces and 
torques that act along the filament when the filament ro- 
tates. For large friction of the load particle, when its size 
is comparable to a bacterial cell body, buckling is mostly 
due to the thrust force created along the filament and 
similar to a rod that buckles under its own weight. In 
the limit of small friction of the load particle, the critical 
thrust force tends to zero and buckling is mostly driven 
by the frictional torque acting along the filament. How- 
ever, our modeling reveals that subtle details of the spe- 
cific problem are important. One has to take into account 
the precession of the buckled filament about the applied 
torque and, in particular, a perpendicular frictional force 
due to the motion of the model bacterium that stabilizes 
the filament against buckling. Finally, taking into account 
the small deformation of the rotating helical filament, we 
are able to obtain a quantitative agreement with the simu- 
lated graph, critical force versus torque, in the biologically 
relevant regime. 

To further illustrate the biological relevance of the ob- 
served buckling transition, we first summarize a few exper- 
imental values. Hotani gives the motor torque for observ- 
ing a polymorphic transition of the flagellum at around 
l.lpNfjm [14], whereas Darnton et al. mention a mean 
torque acting on a flagellum of about 1.4pN|am [5]. These 
values agree with the torques where we observe buckling 
for realistic cell body sizes (see fig. [7]). Reference [5] also 
mentions the relative stiffness of the helical filament so 
that it hardly deforms under rotation which agrees with 
our simulations. Finally, thrust forces created by the bun- 
dle are given as 0.41±0.23pN [5] or 0.5pN [29]. This agrees 
with an estimate F = Gnrjav w 0.6pN where we take the 
radius of the load particle as a = l\im and use the swim- 
ming velocity v = 30|Lim/s. All these values are close but 
below the simulated values F& ~ lpN for real cell sizes. 
However, we note that F c \ scales as A e ^/L^ as our an- 
alytic model shows, and thereby depends on the explicit 
choice of the bending (A) and torsional (C) rigidities. We 
have chosen particular values for them and also set A = C 
in our simulations. So F c \ will vary with the actual param- 
eters. 

It is clear that swimming bacteria should avoid buck- 
ling for efficient locomotion. However, they cannot simply 
increase bending rigidity A since a certain flexibility is 
necessary during polymorphic transformations or when a 
bundle forms. Reference [6] shows pictures where single 
flagella are in a bent conformation similar to the buckled 
state in our simulations. This might be a hint that flagella 
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naturally buckle under their own thrust. In peritrichous 
bacteria such as E.Coli and Salmonella, several flagella 
form a bundle which then has larger bending stiffness and 
therefore buckling is not observed. 

Monotrichous bacteria only use a single flagellum. Their 
conformation differs in pitch and radius from the flagella of 
peritrichous bacteria [48]. A detailed analysis shows that 
their swimming efficiency is reduced due to a smaller pitch 
angle with sin a w 0.75 [49]. This increases the critical 
force F cr oc 1/sina by about 10% compared to peritric- 
hous bacteria and might be an adaption of the monotric- 
hous bacteria to enhance the stability of their single flag- 
ellum. 

We also showed that a pulling flagellum is not sta- 
bly aligned along the applied torque. So most bacteria 
use their flagella to push themselves through the fluid. 
Nevertheless, there are some marine bacteria that use a 
back- and- forth rather than a run-and-tumble strategy for 
chemotaxis. They live in a turbulent aqueous environment 
in the ocean where they experience large shear gradients 
on the micron scale [50]. Simulations in Ref. [50] show 
that in addition to the shear-driven reorientation of the 
bacterium there must be further contributions to the re- 
orientation. Besides rotational diffusion this could also be 
the unstable orientation of the rotating filament when it 
pulls the cell body. Recent experiments on the back-and- 
forth motion of marine bacteria Vibrio alginolyticus di- 
rectly show this reorientation of the flagellum [51]. 

We close with this comment and hope that our work 
initiates a more careful search for the buckling transition 
in bacterial flagella. 

We thank M. Graham and R. Netz for stimulating discus- 
sions and acknowledge financial support from the VW foun- 
dation within the program "Computational Soft Matter and 
Biophysics" (grant no. 1/83 942). 



A Summary of resistive force theory for a 
helix 

At low Reynolds number the force F and torque M acting 
on a particle of arbitrary shape are linearly related to its 
translational and rotational velocities [52] , 



F 

M 



A C 

C T B 



(32) 



The translational friction tensor A, the rotational friction 
tensor B, and the coupling tensor C are determined by 
the shape of the particle. Note that the rotational fric- 
tion tensor B and the coupling tensor C depend on the 
choice of the origin of the coordinate system whereas the 
translational friction tensor A is unique. 

In a moving helical filament, different parts interact via 
hydrodynamic interactions. Nevertheless, using slender- 
body theory, Lighthill demonstrated that one can describe 
the hydrodynamic friction of the filament with the help of 
resistive force theory [40 ,41 . In this theory one introduces 



local friction coefficients per unit length parallel and 
perpendicular (7^) to the tangent vector of the filament. 
Lighthill adjusted the coefficients for the helical filament 
to [40] 



7|| 



27TTJ 



ln(2 ? /r) 



and 7j_ 



In(2g/r) + l/2" 



(33) 



Here rj is the shear viscosity, r = 0.02 |Ltm the cross-sectional 
radius of the bacterial flagellum, and q a characteristic 
length, for which Lighthill derived q = 0.09vl, where A — 
2ttR/ cos a is the filament length of one helical turn. 

In a helical filament with translational velocity v and 
angular frequency uj each segment moves with a velocity 
v + uj x r, where r is the position vector from a point 
on the helical axis to the segment. The force and torque 
densities to initiate such a motion are 

/= (7|| p || +7±P±) •(Hwxr) (34) 
m = r x f, (35) 

where we use the projectors on the local tangent vector 
e 3 and the space perpendicular to it, 



P11 = e 3 e 3 , and 



e 3 e 3 . 



(36) 



Integrating force and torque densities along the helical 
filament with position vector 



r(s) 



( ( cos a 
( it cos ( — 



R 



\ ^ . /cos a \ . \ T 
j , R sin y— ^— sj , sin as J , (37) 



gives Eq. (32). For comparing theory and simulation in 
sect. [3] we calculated the integrals using the computa- 
tional software program "Mathematica" . In particular, we 
took into account that the helical filament in the simula- 
tions does not consists of an integral number of helical 
turns and that the rotational axis is shifted against the 
helical axis. In our analytical theory for the buckling tran- 
sition in sect. [4] we used friction coefficients calculated for 
a full helical turn with L = 2ttR/ cos a. The relevant co- 
efficients become 



A± =Lj ± [l 



cos 2 a 



A\\ = Lj ± (l + £ sin 2 a) 

B\\ =L<y ± R 2 (1 + ^(1 -sin 2 a)) 

Cm = Lj^RS sin a cos a. 



(38a) 

(38b) 
(38c) 
(38d) 



where we use S 



7||~7-L 
7^ 



to characterize the anisotropy 



in the local friction coefficients. Note that An also holds 
for arbitrary filament lengths when L is not an integer 
of a full helical turn. For all other coefficients one obtains 
corrections of the form sin (^pL) that vanish in the limit 
L — > 00. 

The effective friction coefficients for the helical rod 
used in Sec. [4] follow by dividing the friction coefficients 
of Eqs. ([38]) by the rod length L = Lsma: 



a ± 



A±/L , 
B\\/Lq, 



a \\ 

C|| 



(39a) 
(39b) 
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B Rotational motion of a rigid helix 



Starting from Eq. (32) we set v = and concentrate on 



the rotational motion due to a constant external torque 
M with the relevant equation 



M = Buj. 



(40) 



The rotational friction tensor B is symmetric. In the fol- 
lowing we use its frame of eigenvectors {ei, e 2: es} and 
the eigenvalues Bi, B 2 , and B%. We differentiate Eq. ( |40| ) 
with respect to time t, use dfSi = uj x e^, and obtain in 
the eigenframe of B, 



Bid t uJi = (B 2 - B 3 )uj 2 uj 3 , 
B 2 d t uJ 2 = (B 3 - £1)^1^3, 
B 3 d t uJ 3 = (Bi - B 2 )uj x uj 2 . 



(41) 
(42) 
(43) 



These equations are the same as the Euler equations for 
a rigid body with an inertia tensor B and without fric- 
tion. The external torque is zero so that angular momen- 
tum is conserved. Following this analogy and according to 
Eq. (40), the constant external torque in our case corre- 



sponds to the angular momentum of the rigid body, and 
the dissipated energy P — M ■ uj to the rotational kinetic 
energy. Hence, besides the square of the applied torque 
M 2 = Y2i Bfu 2 also the dissipated energy P = J2i BiUjf 
is a conserved quantity and the trajectory of uj follows 
from the intersection of two ellipsoids as illustrated in fig. 
[lOja). In particular, if two of the friction coefficients B{ 
are equal, the angular velocity uj precesses in real space 
on a cone about the direction of the torque [42|l43] . 

We already calculated one component of the rotational 
friction tensor B of a helix with filament length L in the 
previous section. In general, for a long slender helix like 
the normal form of the bacterial flagellum two eigenvalues 
of B are equal to a good approximation, B\ « B 2 oc L 3 . 
The third small friction coefficient B 3 oc L belongs to the 
principal axis, which is parallel to the helical axis, again 
to a good approximation. Hence, a rigid helical filament 
precesses about the applied torque [fig.flOFb)] and does not 
align parallel to the torque as observed in our simulations. 



C Effective bending rigidity of a helix 

We aim at replacing the helical filament by a rod with 
an effective bending rigidity A e ^. Our strategy is to apply 
a small constant torque perpendicular to the helical axis, 
rewrite the total elastic energy as a function of the torque, 
and compare this result with the case of a simple rod to 
obtain A e R. To bend a simple rod with a constant curva- 
ture i?, one needs the bending energy T = ^A e ^f2 2 Lo. 
Using the torque M — A e ^Q [see, for example, Eq. (20)), 
we obtain 



1M 2 

^ A eff 



(44) 



Now we calculate the corresponding elastic energy for the 
helical filament. We apply a constant torque M = Me x 




(b) 



M 



Fig. 10. (a) The constant applied torque and the dissipated 
energy define two ellipsoids (red and blue) in the body fixed 
frame of a rigid helix. The intersection gives the path of the 
angular velocity uj. For a long slender helix two directions are 
degenerate and the trajectories are circles, (b) In the lab frame 
the helix rotates about its axis which precesses about the ap- 
plied torque M. 



perpendicular to the helical axis e z and replace in Kirch- 
hoff 's energy density Q the components of the angular 
strain vector 17 by the components of the torque M = 
AQiei + A(f2 2 - K)e 2 + C(J? 3 - r)e 3 : 



f -I [Mi 

2\ A 



M 2 



c 



(45) 



Note that the components of the applied torque, Mi = 
Me x - Si, depend on the local material frame {ei^e 2 ^es} 
of the helical filament. In leading order in M, we calculate 
the components Mi for the undeformed helical filament 
of Eq. 37 using the Frenet frame and integrate Eq. rt45h 
along the filament: 



m 2 l n 



sin 2 a 



cos 2 a 
C 



O 



sin 2kL\ 
2kL ) 



(46) 



where we used k = cosa/R. We compare this result with 
Eq. (44) and introduce the helix height Lq = Lsina in 



order to identify the effective bending rigidity 



1 



1 1 1 

2 sin a A 



sin 2 a 



A 

C 



cos 2 a 



(47) 



To verify the applicability of the effective bending rigid- 
ity, we study in detail the reorientation rate A of the 
fixed flagellum, when the thrust force pulls at it (see sect. 
3.1.2). Our claim is that the reorientation rate depends 



on the bending of the helical filament as a whole and thus 
A e ff should be the relevant parameter. We therefore deter- 
mined A as a function of the motor torque M for different 
values of the bending rigidity A and the torsional rigidity 
C (see fig. 11). In addition to the helical geometry of per- 
itrichous bacteria used in this paper, we also considered 
the flagellum of monotrichous bacteria, which has different 
helical parameters: = 2. 2 /urn and tq = 2.5/[xm [48] . 
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= A A/C= 1 


• A 


= 2A A/C= 1 


■ A 


= 2A Q A/C = 2 


A A 


= 0.5 A) A/C = 1 


A 


= \A Q A/C = \ 



(ii) 



t A = A A/C= 1 
o A = A A/C = 2 



M , 



Aeff\ 

L J 



Fig. 11. Reorientation rate A as a function of motor torque 
M for different elastic constants, where Aq — 3.5pNu.m 2 is the 
bending rigidity used in this article. Two helical geometries of 
the flagellum are considered: (i) for peritrichous bacteria used 
in this article and (ii) for monotrichous bacteria. 



Dimensional analysis suggests to rescale the torque M by 
the characteristic bending moment Aqr/Lq as in Sect, g] 
and the reorientation rate A by A e ^/{\L^A_[_) ? wh ere A±_ 
is the friction coefficient introduced in Eq. (38a). With 



such a rescaling all different curves for the reorientation 



rate fall onto a common master curve in fig. 11 The effec- 
tive bending rigidity A e ^ is therefore the right parameter 
in an effective description of the helical filament. 
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